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Abstract 

Optimal controller design is addressed for a class of linear, time-invariant systems which 
are dissipative with respect to a quadratic power function. The system matrices are assumed 
to be affine functions of uncertain parameters confined to a convex polytopic region in the 
parameter space. For such systems, a method is developed for designing a controller which 
is dissipative with respect to a given power function, and is simultaneously optimal in the 
linear-quadratic-Gaussian (LQG) sense. The resulting controller provides robust stability as 
well as optimal performance. Three important special cases, namely, passive, norm-bounded, 
and sector-bounded controllers, which are also LQG-optimal, are presented. The results give 
new methods for robust controller design in the presence of parametric uncertainties. 


1 Introduction and Problem Statement 

Over the past three decades, a number of design methods have been developed for control of 
dynamic systems with uncertainties. The methods include H^, y- synthesis, passivity, etc. 
The objective is to obtain closed-loop robust stability and optimal performance in spite of 
model uncertainties. A large class of dynamic systems can be characterized as “dissipative 
systems” which have the property that some of the energy put into the system gets dissipated. 
The term “energy” is defined in a very general sense, and thus the property of dissipativity 
^Senior Scientist, Dynamics and Control Branch 
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encompasses a number of subclasses (for example, passive systems). A major advantage of 
dissipative systems is that they can be robustly stabilized by a controller that itself satisfies 
a certain dissipativity condition. 

This paper presents an approach for the control of dissipative systems via linear quadratic 
Gaussian (LQG) controllers that are also restricted to be dissipative. The method employs 
the above-mentioned stabilization result to ensure robust closed-loop stability, and optimal 
LQG-type controller to achieve the required performance. The approach generalizes the 
LQG-optimal passive controller design of [1] to a broad class of dissipative systems. The 
method is subsequently specialized to three important subclasses of dissipative systems, 
namely, passive systems, norm-bounded systems, and sector-bounded systems. 

1.1 Dissipative Systems 

Consider a linear, time-invariant system E: 


x = Ax + Bu ; y = Cx 


( 1 ) 


where x, u, y are n-, m-, and m- dimensional state, input, and output vectors respectively. 
We assume that A has all eigenvalues in the closed left-half plane. 

Define the quadratic “power function” 


p(y , «) = [y T ’ yT . 


Q N 


y 

iW R 


u 


The following definition is based on [2], [3]. 


(2) 


Definition- The system E is said to be dissipative with respect to the power function p(y, u ) 
if there exists a positive definite quadratic “energy function” E(x) = x T Px with P = P T > 0 
such that the following dissipation inequality is satisfied: 

[ p(y,u)dt>E[x(T)]-E[x( 0)] (3) 

Jo 


VT e [0, oo) and G C™ e . 
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(£!>e denotes the extended Lebesgue space of functions that are square-integrable over all 
finite intervals). The dissipativity condition above can also be expressed in the differential 
form as: 

J t E (x) <p(y,u) (4) 

Three important special cases of dissipative systems are defined below. 

Definition- A system which is dissipative with respect to the power function p(y, u ) in Eq. 
(2) is said to be 

• Passive if Q = 0, R = 0, and N = I. 

• Norm-bounded if Q = —I, R = 7 2 J, and N = 0 for some finite 7 > 0. In this case, 
7 > ffoo-norm of the system. 

• sector-bounded inside the sector [0,6], a < 0 < b if Q = —I, R = — abl , and N = al 
with a = (a + b)/2. 

It should be noted that norm-bounded and sector-bounded systems have to be stable (i.e., 
all eigenvalues of A must have negative real parts) while passive systems can be marginally 
stable. 

1.2 Stability of Feedback Interconnection 

The following result, referred to as the dissipativity lemma, is from [4], 

Lemma- The system E is dissipative with respect to the power function p(y, u) in Eq. (2) 
iff there exists a symmetric, positive-definite matrix P and matrices L and W such that the 
following equations are satisfied 

A T P + PA = C t QC - L t L 
PB = C t (QD + N) - L t W 
R + N t D + D t N + D t QD = W T W 
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(5) 

(6) 
(7) 



Or equivalently, iff the following linear matrix inequality (LMI) is satisfied for some sym- 
metric positive definite matrix P 


A T P + PA- CFQC PB - C t (QD + N) 

B T P - ( QD + N) t C ~(R + N t D + D t N + D t QD) 


(8) 


A definition of strict dissipativity is given next, which is a generalization of the “weak 
strict positive real (WSPR)” definition given in [5]. 


Definition- If A is Hurwitz, (A, L) is observable, and [A,B,L,W] is minimum-phase, the 
system E is said to be strictly dissipative with respect to power function p(y,u). 

We consider feedback-interconnected dissipative systems Ej. and S 2 , which are dissipative 
with respect to power functions pi and p 2 defined as follows: 


Pi(y,u) = [y T ,u T [ 


1 

<0 

£ 

1 


y 

ft? 


u 


The following result from [4] gives a sufficient condition for stability. 


(9) 


Theorem 1- Consider two systems Ei and S 2 connected in the negative feedback config- 
uration (Figure 1). Suppose Ei is dissipative with respect to pi and S 2 is dissipative with 
respect to p 2 . Then the interconnected system is Lyapunov stable if there exist scalars oq > 0 
and a 2 > 0 such that 


onpi{y,u) + a 2 p 2 (-u,y ) <0 Vu,y G R m (10) 

Furthermore, if at least one of the systems is strictly dissipative, then the interconnected 
system is asymptotically stable. 

The systems Ex and S 2 can be considered to be the plant and the controller, respectively. 
The significance of above result is that the stability is robust to model uncertainties; as long 
as the plant is dissipative with respect to power function p \ , any controller which is dissipative 
with respect to power function p 2 will stabilize it, provided that the condition of Theorem 1 is 
satisfied. For example, if the plant is passive, it is stabilized by any strictly passive controller 
[5]; if the plant is norm-bounded by 7, any controller which is strictly norm-bounded by I/7 
will stabilize it; and if the plant belongs to sector [a, b ], any controller which belongs (strictly) 
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to sector — ^] will stabilize it. In the case of passive systems, it is possible to further 
weaken the requirement of strict passivity to marginally strict passivity [6]. 

We consider a class of systems S(9) given by 

x = A{9)x + B(9)u + v; y = C{9)x + w (11) 

where 9 E R k denotes the vector of uncertain parameters and x, u, y are n-, m-, and m- di- 
mensional state, input, and output vectors respectively, v and w denote zero-mean Gaussian 
white noise processes with covariance intensities Qf and Rf respectively. A(9), B(9), C(9) 
are appropriately dimensioned matrices that are assumed to be affine functions of the pa- 
rameter vector 9. The parameter vector 9 assumed to lie in a convex poly topic region V in 
the parameter space, bounded by vertices , j = 1, 2, • • • l. For example, if each component 
of 9 lies in an interval V would be a hyper- rectangular box in R k , and the number 

of vertices l = 2 k . We assume that the system is stable (i.e., A(9) has all eigenvalues in the 
open left-half plane), and the realization is minimal in the entire region V . The system is 
assumed to be strictly proper because the objective is to design an LQG-type controller. 

Suppose the nominal system is represented by the following minimal realization 

x = Ax + Bu + v ; y — Cx + w (12) 


The problem is to obtain a controller which minimizes the performance function 

1 r T 

J = lim — £ / (. x T Q r x + u T R r u)dt (13) 

T->oo ± Jo 


and maintains closed-loop stability in the presence of parametric uncertainties, where Q r > 0 
and R r > 0 are appropriately dimensioned symmetric matrices. 


1.3 Dissipativity in the Presence of Uncertain Parameters 

The following theorem gives a sufficient condition for the system to be dissipative in the 
presence of parametric uncertainties, and is a generalization of the result in [7] for passive 
systems. 
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Theorem 2-The system [A(0), B(9), C(0)] is dissipative with respect to the quadratic power 
function p(y, u ) defined in Eq. (2) V9 e V if there exists a matrix P = P T > 0 such that 


z(e j ,p ) 


A(9 j ) T P + PA(9 j )-C(9 j ) T QC(9 j ) PB(9 j ) - C(9 j ) T N) 

B(9 j ) T P — N T C(9 j ) -R 

for j = 1,2, ---l (14) 


Proof- The proof is similar to that for the quadratic stability case [8]. Because V is a convex 
region in the parameter space, the function P = x T Z(9,P)x is convex in 9 \/x £ R n+m . 
Therefore P takes on its maximum value at one or more of the vertices of V, and thus 
Z{9, P) < 0 M9 £ V if (14) is satisfied. ■ 

The significance of this property (“robust” dissipativity) is that any controller satisfying 
condition (10) can robustly stabilize the system. For example, for passive systems, one can 
determine the largest region V in which the system remains passive [7]. For norm-bounded 
systems, one can determine the largest V in which the system norm remains below a certain 
value (7), or for given V , the smallest norm of the system. For sector-bounded systems, one 
can determine the largest V for a given sector [a, b], or the smallest sector for a given V. For 
these three cases, any controller that is (respectively) passive, or norm-bounded (by 7 -1 ), or 
inside the sector [— b~ l , —a -1 ], will robustly stabilize the system. The next section addresses 
design of such robust controllers that are also optimal. 


2 Optimal Controller 

Given that an uncertain plant is dissipative with respect to a quadratic power function pi, 
the approach considered in this paper is to design an optimal controller which is restricted 
to be dissipative with respect to a quadratic power function p- 2 chosen in such a way that 
the conditions of Theorem 1 are satisfied. 
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2.1 A Special Realization 


Suppose the system (12) is dissipative with respect to the quadratic power function pi with 
weights (Qi,Ni,Ri). Then 3 a matrix P 1 = P[ > 0 such that (8) is satisfied. Suppose 

r e R nxn is a square root of Pi, i.e., Pi = F T r. Using the coordinate transformation 

£ = Tx, the system (omitting the noise terms) becomes 

£ = T£ + Bu (15) 

y = C'£ (16) 

where A = TTr -1 , B = TB, C = CT _1 . Premultiplying (8) (with D = 0) by diag[T ~ T , I m ] 
and postmultiplying by its transpose, we have 

A t + A-c t q 1 c b-c t n 1 ] , = _\q a Qi 2 ] <0 

B T — NjC -Pi ’ Qf 2 Pi 

Therefore, without loss of generality, it will be assumed that the system is in this form, 
i.e., (8) is satisfied with P = I, and the “hat” notation will be dropped in the subsequent 
material. 


2.2 LQG-Optimal Controller 

For the nominal system, it is well-known that the controller which minimizes J consists of a 
linear-quadratic regulator (LQR) and a Kalman-Bucy filter (KBF), and is given by: 


x c = A c x c + B c y (18) 

y c = C c x c u = -y c (19) 

where 

A C = A- BR; l B T P r - P f C T Rj l C B c = P f C T Rf C c = R~ l B T P r (20) 

P r A + A T P r - P r BR~ 1 B T P r + Q r = 0 (21) 

P,A t + AP, - PjC T RpCPf + Q f = 0 (22) 
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The problem is to find an LQG controller that is dissipative with respect to the quadratic 
power function p- 2 with given weights (Q 2 , N 2 , R 2 ). It is assumed that Q - 2 < 0, which is 
usually the case for most dissipative systems. 

Two cases, R 2 > 0 and R 2 = 0, will be considered separately. The following theorem 
gives an LQG controller that is dissipative with respect to p- 2 for the case R 2 > 0. 

Theorem 3- Suppose R 2 > 0, and the LQG performance function weights are such that 
R r 0) Qr ^ d? and 

Rr > N 2 R^ l Nj - Q 2 (23) 


Q r > P r [C T Rj 1 Rj 1 Rj 1 C - C T Rj 1 Rj 1 N^Rj 1 B T - BRj 1 N 2 Rj' Pj'C]R. 

-C T Rj l CP r - P r C T Pj l C 

Q f = Q a + C t {RJ 1 -Ch)C >0 


(24) 

(25) 


Then the resulting LQG controller is strictly dissipative with respect to the power function 


piiih ») = W ,« T 


1 

o 

to 

to 

1 


y 

_ N? Ri _ 


u 


(26) 


Proof- The controller (T c , B c , C c ) is dissipative with respect to p 2 iff 3 a matrix P c = Pj > 0 


such that 


A T C P C + P C A C - CjQ 2 C c P c B c - CjN 2 


BjP c - NjC c 


< 0 


-Pa 


Since R 2 > 0, the above is equivalent to the following condition [9] : 


(27) 


AlP, + P,.A, -CfQ 2 C c +(P c B c - CjN,)RA i (B'P c - NjC c ) < 0 (28) 

After substituting for A c , B c , C c , (28) becomes 

(A - BR~ 1 B T P r - PfC T Rj l C) T P c + P C (A - BR; l B T P r - P f C T R~ f l C) 
-P r BR; 1 Q 2 R; 1 B T P r 

-{PcPfC^R-f 1 - P r B R ~ 1 N 2 ) f?2 1 (P c Pf C T PJ 1 - P r BR~ l N 2 ) T 

< 0 (29) 
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Noting from (17) and (25) that the filter algebraic Riccati equation (ARE) (22) is satisfied 
with Pf = I, using the control ARE (21), and letting P c = P r , (29) can be written as 


—Q r - P, B[R;' + r;'Q 2 r- ] - r,;'n 2 rj'nJ R;']B T p r 
+P r [C T R- f l R^ l R- f l C - C T Rj'R 2 'Nj R; l B T - BR; 1 N 2 R^ 1 Rj 1 C]Pr 

— C T Rj 1 C P r - P r C T Rj l C < 0 (30) 


From (23), 


R r + eh ~ N 2 R^ l Nj > 0 


Pre- and post-multiplying the above by R r L , 


r; 1 + r; 1 q 2 r; 1 - r; 1 n 2 r 2 1 n^ r; 1 > o 


(31) 


(32) 


From (24) and (32), (30) [and therefore (28)] is satisfied (with strict inequality “< ”) and 
the LQG controller (A c ,B c ,C c ) is dissipative with respect to p 2 (y,u). 

Inequality (27) holds in the strict sense; hence in (27), 


Z c =- 



[L c W c ] < 0 


(33) 


where L c e R p/n and W c € R p/m (p > n + m) are of rank n and m respectively. Because 
(27) holds in the strict sense and Q 2 < 0, A c is Hurwitz. ( A c , L c ) is observable because L c 
is of rank n. The transfer function of ( A c , B c , L c , W c ) is non-square and has no transmission 
zeros. Thus the controller is strictly dissipative with respect to the power function p 2 . ■ 

Remark 1- For most common dissipative systems, usually Q i < 0; therefore, (25) is 
satisfied for any Rf >0. 

Remark 2- It should be noted that it is always possible to choose the performance 
function weights to satisfy the conditions of Theorem 3. The procedure for designing an 
LQG controller that is dissipative with respect to a given power function p 2 consists of the 
following steps: 


• Obtain a solution Pi to the dissipativity LMI for the nominal system and transform 
to the special realization of (15), (16) 
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• Choose LQR control weighting matrix R r > 0 satisfying (23), and the state weighting 
matrix Q r based on the performance requirements 


• Solve the control ARE (21) to get P r 

• Choose Rf such that (24) is satisfied; choose Qf as in (25) 

• The controller is given by (20) with Pf = I 


Remark 3- It should be noted that the KBF weighting matrices Q f and Rf are used as 
design parameters and have no statistical significance. 

When R -2 = 0, perhaps the only case of interest is when Q 2 = 0 and N 2 = /, which 
corresponds to passivity. In this case, the LQG controller is passive iff 


A T C P C + P C A C P C B C — Cj <q 

BjP c -C c 0 J ~ 

As shown in [1], the above condition is satisfied with P c = P r if 


(34) 


Q r > BR r l B T 


(35) 


Rj = R, (36) 

Q f = ~{A + A T ) + BR~ l B T (37) 


It should be noted that the above LQG controller is weakly strict positive real (WSPR), 
i.e., it is stable, minimum-phase, and positive real, even when the open loop system is only 
marginally stable. 

When R -2 >0, two important special cases are norm-bounded and sector-bounded LQG 
controllers. 

Corollary 3.1: Norm-bounded controller- Suppose the LQG performance functions 
weights are such that Q r > 0 and 

Rr > I 

Q r > -,- J P t C T B f -CP, - C'Rf'CP, - P/fR/C 
Q, = -iA + A' \+C r ll,'C 
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(38) 

(39) 

(40) 



Then the controller is norm bounded by 7 , i.e., ]|G' c (s )|| 00 < 7 , where G c (s) = C c (sl — 
A c )~ l B c . 

Note that (38), (39), and (40) correspond to (23), (24), and (25) of Theorem 3 (with 
Q-2 = -I, R-2 = 7 2 /, N 2 = Q). 

Corollary 3.2: Sector-bounded controller- Suppose the LQG performance functions 
weights are such that Q r > 0 and 


R r > 



(fl + fe) 2 ' 

4ab 


I 


(41) 


Qr>- Q) P r [C T RfC - {C t RJ 1 R; 1 B t + BR~ 1 Rj 1 C)(a + b)/2]P r 

-C T R~ f l CP r - P r C T Rj l C (42) 

Q f = ~{A + A t ) + C T R~ f l C (43) 

Then the controller belongs to the sector [a, b] (a < 0 < b) in the strict sense. 

Remark 4- Design of LQG controllers that belong to the [0,-1) sector for systems 
belonging to the [a, 00 ) sector was considered in [ 10 ], [ 11 ]. 


3 Concluding Remarks 

Design of linear quadratic Gaussian (LQG) controllers was considered for a class of uncer- 
tain systems which are dissipative with respect to a quadratic power function. The system 
matrices were assumed to be affine functions of parameters belonging to a convex polytopic 
region. A method was given for designing an LQG controller that is dissipative with respect 
to a given quadratic power function. Three important special cases, passive, norm-bounded, 
and sector-bounded controllers, were presented. By appropriately choosing the weighting 
functions, the controller can be designed to provide optimal performance as well as robust 
stability in the presence of parametric uncertainties. 
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Figure 1: Interconnected dissipative systems 
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